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Abstract 

Tumor growth has a number of features in common with a physical process known as molecular 
beam epitaxy. Both growth processes are characterized by the constraint of growth development 
to the body border, and surface diffusion of cells/particles at the growing edge. However, tumor 
growth implies an approximate spherical symmetry that makes necessary a geometrical treatment 
of the growth equations. The basic model was introduced in a former article [C. Escudero, Phys. 
Rev. E 73, 020902(R) (2006)], and in the present work we extend our analysis and try to shed light 
on the possible geometrical principles that drive tumor growth. We present two-dimensional models 
that reproduce the experimental observations, and analyse the unexplored three-dimensional case, 
for which new conclusions on tumor growth are derived. 

PACS numbers: 87.10.+e, 68.35.Fx, 87.19.-j 
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I. INTRODUCTION 



One of the highest mortahty rates in developed countries is due to cancer. For this reason, 
it is one of the most studied diseases, and however, it is still far from being well understood. 
Aside from the intensive medical and biological research, an increasing number of theoretical 
models is being introduced in order to describe some of the fundamental properties of tumors. 
One of the most common mathematical approaches is the use of partial differential equations, 
which has lead to some interesting results in the field l|. 

A different methodology was used by Bru et al. 0, 3|, that employed some tools of fractal 
geometry, as scaling analysis, to characterize the rough interface of growing solid tumors. 
They found strong empirical evidence that a broad class of tumors belong to the same 
universality class: the molecular beam epitaxy (MBE) universality class. MBE is a well 
known process in physics, in which a crystal surface grows due to the external input of 
atoms coming from a beam directed to the growing surface . This finding allowed them to 
postulate that tumor dynamics has some features in common with MBE, say: (a) a linear 
growth rate, (b) the constraint of growth activity to the outer border of the tumor/crystal, 
and (c) surface diffusion at the growing edge. All of these features were again tested against 
experiment with a positive result. 

These experimental observations lead to a new and very interesting picture of tumor 
growth. The usual assumed exponential growth is replaced by a linear growth, and the 
dynamics is constrained to the peripheral region, because it is assumed that what happens 
in the core of the tumor has little effect on growth. Surface diffusion has been identified as 
a mechanism for favoring tumor growth. The host tissue exerts pressure on solid tumors 
which opposes their growth, but surface diffusion drives the tumor cells to the concavities 
of the interface, keeping the number of neighboring cells that belong to the host tissue 
to a minimum. Since these cells are responsible of the pressure exerted on the tumor, 
surface diffusion minimizes the pressure on the interface, and favors the propagation of the 
tumor 

This new theoretical description of a solid tumor was used to develop a strategy to stop 
the growth [sj. Suggested by the above findings, Brii et al. proposed that an enhancement of 
the pressure on the tumor surface would be able to decelerate, and eventually to stop tumor 
growth. They performed a new experiment in which they observed the response of the tumor 
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to an enhancement of the inmune response. An increase of the number of neutrophils shifted 
the dynamics of the interface from the MBE universahty class to the much slower quenched 
Edwards- Wilkinson (QEW) universality class, and to the pinning of the tumor interface P|. 
This technique was later applied to a patient with a terminal cancer, who subsequentialfy 
improve his state till finally achieve a good health, possibly due to the applied treatment p . 

This success has a fundamental importance for the development of efficient therapies, 
and therefore it would be highly desirable to achieve a good theoretical understanding of 
the models used. MBE dynamics is described by the MuUins-Herring equation 0, 0] 

dth = -KV^h + F + r]{^,t), (1) 

where h is the interface height, K is the surface diffusion coefficient, and ri{'K, t) is a Gaussian 
noise with zero mean and correlations given by 

< ?7(x, t)r7(x', f) >= D(5(x - x.')6{t - t'). (2) 

This equation was developed to describe crystal growth, so it assumes that the substrate is 
planar and does not change its size in time. This no longer applies to tumors, since they 
are approximately spherical and grow linearly in time; both features are possibly the main 

n 

discrepancies between standard MBE dynamics and the dynamics of tumors [3]. The same 
criticism applies to the equation describing QEW dynamics 0| 

dth = KV^h + F + r]{x,h), (3) 

where ri{'x, h) is a quenched disorder with zero mean and correlations given by 

< r7(x, /i)77(x', h') >= D(5(x - x')A(/i - h'), (4) 

and the function A characterizes the nature of the quenched disorder. 

In a former article, we developed a stochastic partial differential equation describing the 
same dynamics as Eq.(0) but with the correct geometrical symmetries j^. The analysis of 
this equation revealed that it was able to reproduce some of the fundamental mechanisms 
of tumor growth found in experiments . In the present work we extend the geometrical 
approach to tumor growth by developing and analyzing spherically symmetric equations 
describing both MBE and QEW dynamics. T 
experimentally in the (1 + 1)— dimensional case 



lese equations model the behaviour found 
3, 0, 1^ , and allow us to predict what would 



happen in the more reahstic and unexplored case of (2 + 1)— dimensional geometry. We also 
derive the equations using a more systematic technique, that allows us to conjecture what 
are the geometrical principles that drive tumor growth. 



II. EXPANSION FROM A POTENTIAL 

In general, a planar stochastic growth equation may be written in the form 

^^^ = G[h{x,t)]+r^{x,t), (5) 

where G is the deterministic growth term and rj is the noise. If we want the mean field equa- 
tion to describe a conservative dynamics (as for instance a diffusion), then the deterministic 
part must have the form of a continuity equation 

^ = -V-j(..0. (6) 

where the macroscopic current j{x,t) describes the flux of cells on the surface. The current 
j{x, t) arises in general from differences in the local pressure p{x, t) as argued before, following 
the law 

= vn(x,t), (7) 

where 11 is defined as a pressure potential. We can perform the expansion of 11 in terms of 
the pressure 

n(x, t) = -Aip{x, t) + ^2 V^p(x, t) + ■ ■ ■ . (8) 

Since the difference in pressure comes mainly from the differences in height of the different 
parts of the interface, we may assume that p oc h. Finally, we can write the deterministic 
part of the evolution equation as 

= A^V'h - A.V'h + ..., (9) 

where we can identify the terms present in the drifts of both Eq.(^ and Eq.®. 
The equation of growth of a general Riemannian surface reads 

dtf{s, t) = n(s, t)r[r(s, t)] + $(s, t), (10) 

where the d+1 dimensional surface vector r(s, t) = {^^(s, t)}^jl\ runs over the surface as 
s = {s^}f^i varies in a parameter space (in the following, latin indices vary from 1 to d and 



greek indices from 1 to c? + 1). In this equation n stands for the unitary vector normal at 
the surface at r, F contains a deterministic growth mechanism that causes growth along the 
normal n to the surface, and $ is a random force acting on the surface. In our case the 
deterministic part should include a term modelling cell diffusion in the tumor border. When 
surface diffusion occurs to minimize the surface area the corresponding term in the equation 



is Q: 



r, = -kAblH, (11) 



where Abl is the Beltrami-Laplace operator 

Abl = -^d.i^g'^d,), (12) 

Qij is the metric tensor and g is its determinant, di = d/ds^ is a covariant derivative, and 
H = n ■ A^^r is the mean curvature. Summation over repeated indices is always assumed 
along this work. Finally, the unitary normal vector is given by n = g~^^'^difx ■ ■ ■ x ddr. The 
equation of growth can be derived straightforwardly from here [^, but since our aim is to 
understand the geometrical principles underlying tumor growth, we write the more general 
growth equation for those cases in which the drift can be derived from a potential: 

g,rM^.^mhM, (13) 

In our case the potential V depends on the mean curvature H of the interface. This depen- 
dence can be expressed in a power series expansion: 

„ AT N 

From here we can derive straightforwardly the stationary probability distribution functional 
P[r{s, t)], which yields the probability of the interface configuration r(s, t) in the limit t oo 



P[r{s,t)] = A/'exp 



V[r{s,t)] 



(15) 



D/2 

where D is the noise strength and A/" is the normalization constant. In this work, we are 
more interested in the dynamical rather than in the stationary properties of the model, so 
we will focus on the stochastic partial differential equation which describes tumor growth. 
In the general equation, the contribution to the drift reads 

r, = -^h ■^-^ = kJ W^' - iAblH^-' - iw-' V An , (16) 



iAblH^-'~iW-'Y.\^ 
i=i / 



where Xj are the eigenvalues of the matrix of the coefficients of the second fundamental form 
and express the principal curvatures of the surface [l^ . 

III. STOCHASTIC EQUATIONS FOR TUMOR GROWTH 

In this section we will build the models for describing the growth of a non-treated tumor. 
First of all, we will derive the zeroth, first, and second order terms in the expansion Eq. (jl4p 
in one and two dimensions. The corresponding contributions to the drift in polar coordinates 
read 

ro(^=l) = ^|^, (17a) 

Ti{d=l) = 0, (17b) 
1 d^r 

r2(^=l) = --,g^,, (17c) 

1 / d'^r 1 d^r \ 

r.(<' = 2) = ;.(a^ + -i(5)a^j. (™) 

where we have linearized the different derivatives of r about zero and chosen only the most 
relevant terms in the renormalization group sense. Retaining only the linear terms is a valid 
approximation whenever sharp changes in the tumor interface are absent It is very 
important to realize that Ti{d = 1) vanish identically, since H = Xi in d = 1. This is a 
consequence of the Gauss-Bonnet theorem, which states that the integral of the Gaussian 
curvature on a closed surface is a constant. Since H = K in d = 1, the variation of Vi is 
zero. Another important fact is that Ti{d > 1) is strictly nonlinear in the derivatives of r 
in the case of a planar geometry, so it will never survive after a linearization We will 
show what are the consequences of this in the next sections. 

In the case of non-treated tumors, only the r2 term seems to appear in the dynamics. In 
this case, we can derive the following equations describing the tumor interface dynamics in 
(1 + 1)— dimensions 
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where the noise ?7(5, t) is Gaussian, with zero mean, and correlation given by 

< r]{e, t)r]{e', t') >= D6{0 - e')5{t - t'), (19) 

and (2 + 1)— dimensions 

dr__K_ /9V 2 d^r 1 d^r \ 1 

di~ \W + sin2(^) de^dcj)^ + ^ r,/\M^ 



where the noise ri{6, 0, t) is Gaussian, with zero mean, and correlation given by 

< 7]{e, 0, t)7]{e', 0', t') >= D5{e - e')5{(p - (p')6{t - 1'), (21) 

and the noise must be interpreted in the Ito sense. A more detailed derivation of these 

equations can be found in Ref. j^. ^ ^ 

To analyze these equations we will perform a small noise expansion [l^j, where the solution 
is decomposed as follows: 

r{e,t) = R{t) + VDp{9,t), (22) 

in the (1 + 1)— dimensional case, where R is the deterministic solution, given by R{t) = 
Ft + Rq. The stochastic perturbation obeys the equation 

9p K d'^p 1 , , , . 

m = ~{Ro + FtyW + 7^+^"^^ ' ^ ^ 

where the noise Ti[9,t) is Gaussian, with zero mean, and correlation given by 

< r]{e, t)7]{e', t') >= s{e - e')6{t - 1'). (24) 

The discrete Fourier transformed version of this equation reads 

dpn Kn* 1 , . 

where the noise r]n{t) is Gaussian, with zero mean, and correlation given by 

< Vn{t)Vmit') >= i2ny^5n,-mSit - t'), (26) 

where 5n-m denotes the Kronecker symbol. The mean value of the stochastic process obeys 
the equation 

d< Pn> Kn^ ^ ^ 
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that can be solved to yield 



< Pn{t) >= exp 



< Pn(0) > . 



(28) 



We can also calculate the two-point correlation function using standard techniques j3] 

d < Pn(t)pm(t) > = < Pn(t + dt)pm(t + dt) > - < Pn{t) Pm{t) > = 
< dpn{t)pm{t) > + < pn(t)dpm(t) > + < dpn(t)dpm(t) > = 

~ iRo + FtY ^''' + ""'^ ^ Pnit)pm{t) >dt + —l—5^^_^{2n)-'dt, (29) 

where we have used the facts that r]n{t) = {2n)^^/'^dWn{t) / dt and dWn{ty = dt, and where 
dWn{t) denotes the increment of a Wiener process. And thus, the two-point correlation 
function Cn^mit) =< PnifyPmif) > obcys the equation 

dCn,m{t) K 



dt {Ro + Fty 

The exact solution of this equation is 



(n^ + m^)Cn,m{t) + i^^,^ r^, K-m- 



Ro + Ft 



(30) 



Cn,m{t) — GXp 



3F 



R^ (i?o + Ftf 



X 



6nF 



exp 



3FRI 



Ei 



3FRI 



Ei 



'3F{Ro + Fty 



(31) 



where Ei denotes the exponential integral. 

In the (2 + 1)— dimensional case we can perform the expansion 



r{e,(i),t) = R{t) + VDp{e,(i),t), 

where R{t) = Rq + Ft. The stochastic perturbation obeys the equation 



dp 

at 



K 



d^p 



+ 



{Ra + Ftf \de^ sm\e)de 



d^p 1 



d^p 



+ 



1 



{Ro + Ft)^/\M&)\ 



which implies 

dpn,r, 



dt 



K ( A ^ 2 2 ^ a\ ^ / \ 

{Ro + FtY\ 3 3 J^' /?o + ^r 



(32) 



(33) 



(34) 



where A = 160F(7r/4|2)/(637r), and F denotes an incomplete elliptic integral of the first 
kind The noise is again Gaussian, with zero mean, and correlation given by 



< r]n,mit)r]p,q{t') >= {2n) "^dn-pSm-qSit - t'). 



(35) 



We can derive the equation for the first moment 
d < pn,m > _ K 



dt 

and solve it to obtain 



< Pn,m{t) >= exp 



{Ro + Fty 



+ -n^m^ + -m* ) < pn,m >, 



(36) 



K{8m'^ + 8m^n^ + 3n 
9F 



1 1 

{Ro + Ftf\ 



< Pn,m(0) > . (37) 



The equation for the two-point correlation function can be straightforwardly derived and is 



Jt 



K 



{Ro + Fty 



+ 



B 



o^m,—qj (38) 



(i?o + Fi)2""'-^"™'-^' 

where B = 6400F(7r/4|2)V(39 697r4)(Ri 0.03), and Cn,m,p,q{t) =< Pn,m(^)Pp,,(^) >■ The 
solution to this equation is 



Cn,m,p,qit^ — CXp 



X 



-^{^m^ + 8m^n^ + 3n^ + 3p^ + Sp^g^ ^ §^4^) | _i_ 

' 4 , o_2„2 , o„4 , o„4 , o„2„2 , o„4^ 



C'n,m,p,g(0) + i^^ji^—p^^^—q CXp 



/' 



exp 



9FRI 

,2^2 



Bm* + 8m^n^ + 3n^ + 3p^ + Sp^g^ + 



X 



K{%m^ + Sm^n^ + 3n^ + 3p^ + Sp^?^ + 8g 



(Fr + Ra)-^dT 



(.39) 



9F(Ft + i?o)=^ 

One can see that our exact solutions reveal some interesting characteristics of the growth. 
The mean value of the small perturbation decreases in time proving the stability of the mean 
field radially symmetric solution, both in one and two dimensions. The correlation functions 
give a much more interesting information. The correlations generated by the noise are much 
smaller in the two-dimensional case: not only the numerical prefactor is smaller, but in the 
one-dimensional case these correlations decrease in time as an exponential integral, while in 
the two-dimensional case the decay is much faster. This means that the effect of the noise 
is much stronger in one dimension, while the deterministic dynamics are more robust in two 
dimensions, something that might have serious consequences in tumor therapy, as we will 
show below. 



IV. GROWTH IN A DISORDERED MEDIUM 



In the models commented in the introduction, it seems that the disorder is induced by 
the enlarged amount of neutrophils introduced to stop tumor growth. But actually, the 
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medium in which tumors grow is highly disordered. Tumors are extensively infiltrated by 
immune cells which may constitute as much as one third of its volume. Both the tumor 
phenotype and the tumor environment are very heterogeneous. The former is the result of 
accumulating random mutations, variable enviromental selection forces and perhaps restric- 
tion of proliferate capacity in non-stem cell components of the tumor. In addition, the tumor 
environment is extremely heterogeneous primarily due to disordered angiogenesis and blood 
flow. These facts suggest that one possible explanation for the noise term appearing in the 
equations of the last section comes from the underlying disorder. In fact, if we consider 
Eq.(jni) and the corresponding correlations of the quenched disorder Eq.(jH), we see that for 
large F the function h will behave as /i ~ Ft, and assuming that the function A models 
short-range correlations (as found experimentally Isl), implies A(/i — h') oc 5{t — t'). This 
means that, far from the pinning threshold, the role of the disorder is to induce thermal 
fluctuations in the dynamics as those found in Eq.(0). Physically, this means that a rapidly 
moving interface samples so many values of the disorder in a small time interval, that the 
overall effect is that of a time dependent noise. This implies in turn that the intensity of 
the noise is proportional to the disorder, and thus an enhancement in the inmune response 
corresponds to a stronger noise in Eqs. ljlSll^Uj) . while the tumor is still in a state far enough 
from the pinning threshold. 

Not only the noise, but also the diffusion terms vary from Eq.(^ to Eq.Q. We can 
understand this effect if we suppose the terms in the expansion Eq. ()14jl dependent on the 
disorder. Assuming that Kq ~ D and that K2 is independent of where D is the intensity 
of the disorder, then we find that the Kq term is the most relevant in high disorder, while 
K2 would be the most important term in other case. We can use these facts to build a more 
general model of tumor growth in the (1 + 1)— dimensional setting 

To d^r K2 



dr Kn d'^r Ko d^r _ 1 



at r^oe. .4a(,4 + ^' + 7p'K''.''). («) 

where the correlations of the quenched disorder are given by 

< r]{e, r)7]{e\ r') >= 6(9 - e')A{r - r'). (41) 

In weak disorder, this equation reduces to Eq. (|TH|l . because the disorder behaves as a thermal 
noise and the term proportional to Kq loses its importance, as explain above. 

In the (2 + 1)— dimensional case the situation is different. In opposition to the former case, 
now the term proportional to Ki does not vanish identically as shown in section |n] Taking 
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into account this fact, we can again build the more general equation for tumor growth, that 
reads 



dr 
di 



Ko fd'^r 



K2 fd^r 



)02 



+ 



2Ki fd'^r 



+ 



+ F + 



80^ ' sm\e)d(P^ 
:r]{e,(j),r) 



where the correlations of the quenched disorder are given by 

< 7]{e, 0, r)r]{e', 0', r') >= 6(6 - e')6{(f) - (f)')A{r - r'). 



(42) 



(43) 



Now, it is our goal to understand what differences appear in the evolution of the equations 
for tumor growth in different dimensions, mainly due to the presence of the Ki term in the 
(2 + 1)— dimensional model. For this, we will study the deterministic counterparts of Eqs. (|4m 
021) • In the case of (1 + 1)— dimensions we have 

dr Kq d^r K2 d^r 



+ F, 



(44) 



dt r2 Oe^ r4 09^ 

which admits the radially symmetric solution r{t) = Rq + Ft. If we perform the linear 
stability analysis of this solution by adding a small perturbation p, we obtain the equation 



dpn 
dt 

that can be solved to yield 
Pn{t) = exp 



Kon' 



+ 



1 



{Ro + Ftf {Ro + Ft)\ 
Kon^ 



1 



+ 



Pn 



1 



1 



Rq + Ft Ro 



(45) 



Pn{0). (46) 



3F [{Ro + Ft)^ 

This reveals that all the Fourier modes n 7^ of the solution are linearly stable for t > 0. 
The n = mode is marginal, but this is unimportant because it implies perturbations 
homogeneous in 6. On the other hand, the equation corresponding to the (2+1)— dimensional 
case reads 



dr 
di 



Ko fd'^r 



d'^r 



89^ ' sm\e)d(P^ 
/(9V 2 



27^1 /aV 



d'^r 



06^ sm\e)d<i)^ 



f 77777 + 



+ 



36^ sm\e) de^d(f)^ sm\e) 



+ F, 



(47) 



which again admits the solution r{t) = Rq + Ft. If we perform the linear stability analysis 
we arrive at the equation 



dt 



{Ft + Roy 



3 J ' {Ft + Ro)^ 



2 ^2 

n H — m 
3 



K. 



{Ft + RoY 



rY + —rn^ry + -nrY 
o o 



Pn,m 
Pn,my 



(48) 
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that can be solved to yield 

Pn,m{t) = exp 

3F 



K2(8m'^ + 8m^n^ + 3n 



9F 



Rq 



[Rl {Ro + Ft)\ 

Pn,m \}J ) ) 



(49) 



i?o + Ft_ 

and we see again that all the Fourier modes, except m = n = 0, are stable for t > 0. The case 
m = n = is marginal, but it is unimportant since it implies perturbations homogeneous in 
9 and 0. 

The effect of the Ki term can be immediately understood by regarding Eqs. ()4(j^ \^^ . It 
is a new mechanism for "dissipating" curvature. When a stochastic perturbation drives the 
solution away from the radially symmetric form, then this perturbation behaves as stated 
by this two equations. As can be clearly seen, the restoring of the symmetric form is faster 
in the (2 + 1)— dimensional case due to the presence of the Ki term. 



V. CONCLUSIONS 

Motivated by the successful research on tumor growth by Brii et al. , we 

have introduced theoretical models able to reproduce some of the features found in these 
experiments. We have also analyzed these models in order to better understand what is 
happening in the physical phenomenon. The models were derived as an expansion from a 
potential, which resulted a power series expansion in the mean curvature of the surface. The 
effect of the deterministic part of the equations is thus to reduce the mean curvature and its 
different powers, and this way to minimize the pressure in the tumor border and to favour 
tumor growth. 

Since all the experiments were performed with two-dimensional tumors, it is important to 
develop theoretical models that represent what happens in these experiments, and to extend 
them to the three-dimensional case. This way we will get an insight of the evolution of the 
more realistic three-dimensional system, and we can derive conclusions that can be used as 
a guide for future experiments. 

We have analyzed the effect of stochasticity in these equations, and we have shown that 
stochastic effects are much less relevant in the case of (2 + 1)— dimensional growth. Also, we 
have discussed the possible origin of the noise on the underlying disorder of the system, and 
how its intensity can be enhanced by increasing the number of inmune cells in the tumor 
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environment. This shows that the strategy of enhancing the inmune response in order to stop 
tumor growth should be less effective in the case of three-dimensional tumors. In addition 
to this, we have shown that a new term (the Ki term) appeared with higher dimensionality. 
The origin of this term is very interesting, because it vanishes identically in one dimension 
for any geometry and in any dimension in the case of a planar geometry. Thus it appeared 
in the dynamics of the growing tumor as a combined effect of dimensionality and geometry. 
This term contributes to minimize the pressure and to favour tumor growth. Since it is 
present only in the three-dimensional case, it is another new mechanism that helps tumor 
propagation in this dimensionality, what lead us to conclude that it is much more difficult 
to stop a three-dimensional tumor than to stop its two-dimensional counterpart. 

All along this work we have assumed that there are no overhangs in the interface in the 
radial direction, in such a way that a single valued solution of the corresponding stochastic 
growth equation makes sense for representing the interface evolution. There are, however, 
situations for which we cannot assume this. An interesting problem for future work is to 
derive a continuum model allowing overhangs and an arbitrary topology of the growing 
interface. Similar models were developed in a different context [l5[ llfi|. and allow a more 
detailed description of a complex growth phenomenon. 

Acknowledgments 

This work has been partially supported by the Ministerio de Educacion y Ciencia (Spain) 
through Projects No. EX2005-0976 and FIS2005-01729. 



[1] M. Castro, C. Molina-Pans, and T. S. Deisboeck, Phys. Rev. E 72, 041907 (2005). 
[2] A. Bru, J. M. Pastor, I. Fernaud, I. Bru, S. Melle, and C. Berenguer, Phys. Rev. Lett. 81, 
4008 (1998). 

[3] A. Brii, S. Albertos, J. L. Subiza, J. L. Garcia- Asenjo, and I. Brii, Biophys. J. 85, 2948 (2003). 
[4] A.-L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth (Cambridge University 

Press, Cambridge, 1995). 
[5] A. Brii, S. Albertos, J. A. L. Garcia-Asenjo, and I. Brii, Phys. Rev. Lett. 92, 238101 (2004). 
[6] A. Brii, S. Albertos, F. Garci'a-Hoz, and I. Brii, J. Clin. Res. 8, 9 (2005). 

13 



[7] W. W. Mullins, J. Appl. Phys. 28, 333 (1957). 
[8] C. Herring, J. Appl. Phys. 21, 301 (1950). 
[9] C. Escudero, Phys. Rev. E 73, R020902 (2006). 
[10] M. Marsih, A. Maritan, F. Toigo, and J. R. Banavar, Rev. Mod. Phys. 68, 936 (1996). 
[11] Our polar coordinates are, for (1 + 1)— dimensional geometry, r = {r{9)cos{9),r{9)sm.{6)), 
with G [0, 27r]. For (2 + 1)— dimensional geometry they read r = 
(r(6',^)sin(6l)cos(^),r(6',<;^)sin(6')sin(0),r(6',(?!))cos(6')), with ^ G [0,27r] and 9 G [0,7r]. 
[12] J. Garcia-Ojalvo and J. M. Sancho, Noise in Spatially Extended Systems (Springer- Verlag, 

New York, 1999). 
[13] C. R. Doering, Phys. Lett. A 122, 133 (1987). 

[14] The incomplete elliptic integral of the first kind, for — 7r/2 < (j) < tt/2, is F{(l)\m) = Jq[1 — 
msm'^{9)]~^/^d9; see M. Abramowitz and I. A. Stegun (Eds.), Handbook of Mathematical 
Functions with Formulas, Graphs, and Mathematical Tables (9th printing. New York, Dover, 
1972). 

[15] P. Keblinski, A. Maritan, F. Toigo, J. Koplik, and J. R. Banavar, Phys. Rev. E 49, R937 
(1994). 

[16] P. Keblinski, A. Maritan, F. Toigo, R. Messier, and J. R. Banavar, Phys. Rev. E 53, 759 
(1996). 



14 



